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Abstract 

For the spacelike momenta k of the virtual photon 7*, the ^ (p)^/*(k)^/(k') 
transition form factor is considered in the coupled Schwinger-Dyson and 
Bethe-Salpeter approach in conjunction with the generalized impulse approx- 
imation using the dressed quark-photon-quark vertices of the Ball-Chiu and 
Curtis-Pennington type. These form factors are compared with the ones pre- 
dicted by the vector meson dominance, operator product expansion, QCD sum 
rules, and the perturbative QCD for the large spacelike transferred momenta 
k. The most important qualitative feature of the asymptotic behavior, namely 
the 1/k 2 dependence, is in our approach obtained in the model-independent 
way. Again model-independently, our approach reproduces also the Adler- 
Bell-Jackiw anomaly result for the limit of both photons being real. For the 
case of one highly virtual photon, we find in the closed form the asymptotic 
expression which can be easily generalized both to the case of other unfla- 
vored pseudoscalar mesons = 7T°, 77s, 770 > Vc Vb, and to the case of arbitrary 
virtuality of the other photon. Implications thereof for certain important 
theoretical and experimental applications are pointed out. 

PACS: 11. 10. St; 13.40. -f; 14.40.Aq; 14.40.-n; 14.40Gx 
Keywords: Schwinger-Dyson equations; Ball-Chiu and (modified) Curtis-Pennington 
vertices; Chiral symmetry breaking; QCD; QED axial anomaly; Two-photon processes 
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1. The interest in the form factor T 7r o(— Q 2 , 0) for the transition ^f*(k)^f(k') — > 7r°(p) (where 
A; 2 = — Q 2 ^ is the momentum-squared of the spacelike off-shell photon 7*), has again been 
growing lately for both experimental (the new CLEO data |1| and plans for new TJNAF 
measurements [fj]) and theoretical reasons - such as its relevance for the hadronic light- 
by-light scattering contribution to muon g — 2 ||, which is in turn again relevant for the 
experiment E821 at BNL Qj. Nevertheless, this transition form factor still represents a more 
demanding theoretical challenge than usually thought, as the perturbative QCD (pQCD) 
may be not sufficient even at the highest of the presently accessible momenta (Q 2 < 10 
GeV 2 ) - e.g., see Refs. ||. These references therefore indicate desirability of having direct 
calculations of T 7r o(— Q 2 ,0) without any a priori assumptions about the shape of the pion 
distribution amplitude ip n . Even more important for the present work is that the pQCD 
approaches, oriented mostly at reproducing the large Q 2 behavior, have problems at low Q 2 . 
This was most recently stressed in Ref. [[J discussing the status of QCD-based theoretical 
predictions for this process. See especially its Sec. IV for clarifications how such approaches 
can fail to reproduce the Q 2 = value (|1|) corresponding to the Abelian axial alias Adler- 
Bell-Jackiw (ABJ) anomaly which explains the 7r° decay into two real photons (k 2 = k' 2 = 0). 

On the other hand, the approaches (such as the present one) which rely on the chirally 
well-behaved Schwinger-Dyson (SD) and Bethe-Salpeter (BS) equations^ for the light pseu- 
doscalar mesons (71*, K, 77, ...) while taking care that the Ward-Takahashi identities (WTI) of 
QED are respected, satisfy in this respect even what Ref. || calls the "maximalist" require- 
ment - namely that the fundamental axial anomaly result 

should be satisfied automatically in the chiral (and soft) limit. (This is clearly superior 
to imposing it as an external condition, which can by analogy be termed the "minimalist" 
attitude.) In the present approach JTT-^], the Abelian axial anomaly amplitude (|I|) is 



obtained without imposing any requirements or constraints on the solutions for the n wave 



function or the solutions for the quark propagators (e.g., see p4| , p!5| ,|9[|). Since the anomaly 
- and thus also T n o(0, 0) of Eq. (P - must not depend on the internal structure of the pion, 
this is exactly as it should be. What must then be explored in the coupled SD-BS approach 
is if the large Q 2 behavior can be satisfactorily understood in that approach. In particular, 
the comparison with the new CLEO data Q - at Q 2 up to 8 GeV 2 - must be made, since the 
results of the SD-BS approach |Tl[ (as well as the results of the earlier and closely related 
Ansatz approach of [|16[]) have so far been compared only with the older CELLO data [17] 



at Q 2 not exceeding some 2.5 GeV 2 . In this letter we will show how the SD-BS approach 
to modeling QCD provides a good description for both low and high values of Q 2 . In other 
words, it shows how one can in fact obtain something similar to the Brodsky-Lepage (BL) 
interpolation formula T n o(—Q 2 , 0) = (l/47r 2 / 7r ) / (l + Q 2 /8ir 2 f 2 ) proposed |L8j as a desirable 
behavior for the transition form factor because it reduces to the ABJ anomaly amplitude 
(pi) at Q 2 = 0, while agreeing with the following type of leading behavior for large Q 2 : 



1 See, for example, Ref. || ; and the recent reviews ||, or [1C|. 



2 



f 

T^o(-Q 2 , 0) = jj^ (J = constant for large Q 2 ), (2) 



which is favored both experimentally and theoretically (e.g., see the QCD-based predic- 
tions in Refs. IP^8|-|20|,[5|J^,|^] , and references therein.) 

2. In the coupled SD-BS approach, the BS equation for the pion bound-state qq vertex 
T n o(q,p) employs the dressed quark propagator S(k) = \A(k 2 )^ — B(k 2 )]~ 1 , obtained by 
solving its SD equation. Solving the SD and BS equations in a consistent approximation 
is crucial (e.g., see Refs. [pl|-p^|) for obtaining qq bound states which are, in the case of 
light pseudoscalar mesons, simultaneously also the (pseudo-)Goldstone bosons of dynamical 
chiral symmetry breaking (D%SB). 

Following Jain and Munczek |^T],^2j, we adopt the ladder-type approximation sometimes 
called the improved [§] or generalized [[HJ ladder approximation (employing bare quark- 
gluon-quark vertices but dressed propagators). For the gluon propagator we use an effective, 
(partially) modeled one in Landau-gauge |2T| , [2"2"f , given by G(—l 2 )(g PlV — IH V ' /l 2 ) . (This 



Ansatz is often called the "Abelian approximation" f24fl .) What is essential is that the 
effective propagator function G is the sum of the perturbative contribution Guv an d the 
nonperturbative contribution G IR : G(Q 2 ) = G UV (Q 2 ) + G IR (Q 2 ) , (Q 2 = —I 2 ) . The 
perturbative part Guv is required to reproduce correctly the ultraviolet (UV) asymptotic 
behavior that unambiguously follows from QCD in its high-energy regime. Therefore, this 
part must be given - up to the factor 1/Q 2 - by the running coupling constant a s (Q 2 ) 
which is well-known from perturbative QCD, so that Guv is n °t modeled. As in Refs. 
TT|— p~3| , we follow Refs. ||2l|j2^] and employ the two-loop asymptotic expression for a s (Q 2 ). 
For the modeled, IR part of the gluon propagator, we adopt from Ref. |22[] Gir(Q 2 ) = 
(167i 2 /3)aQ 2 e-^ Q2 , with their parameters a = (0.387 GeV)~ 4 and /i = (0.510 GeV)~ 2 . 

More details on our calculational procedures can be found in our Refs. [^T|-|T3[. We 
essentially reproduce Jain and Munczek's [^] solutions for the dressed propagators S(q), 



i.e., the functions A(q 2 ) and B(q 2 ), as well as the solutions for the four functions comprising 
the pion bound-state vertex T^o. Actually, Ref. |22| employs the BS amplitude X-k (1iV) = 



S(q + p/2)T w o(q,p)S(q — p/2), which is completely equivalent. 

3. We assume that the 7r°7*7 transition proceeds through the triangle graph (Fig. 1), and 
that we calculate the pertinent amplitude T^(k,k') = e af3 ' lu k a k' [3 T 7T o(k 2 ,k' 2 ) as in Refs. 
PT^T3[, using the framework advocated by (for example) [H-[T^, 2^|J2^] in the context of 
electromagnetic interactions of BS bound states, and often called the generalized impulse 
approximation (GIA) - e.g., by p6|j25[| . To evaluate the triangle graph, we therefore use 
the dressed quark propagator S(q) and the pseudoscalar BS bound-state vertex Yp(q,p). 
Another ingredient, crucial for GIA's ability to reproduce the correct Abelian anomaly re- 
sult, is employing an appropriately dressed electromagnetic vertex T tl (q',q), which satisfies 
the vector Ward-Takahashi identity (WTI) (q' — q) fl T fi (q' ,q) = S^ 1 (q') — S~ 1 (q) . Namely, 
assuming that photons couple to quarks through the bare vertex 7 M would be inconsistent 
with our quark propagator S(q), which, dynamically dressed through its SD-equation, con- 
tains the momentum-dependent functions A(q 2 ) and B(q 2 ). The bare vertex 7^ obviously 
violates the vector WTI, implying the nonconservation of the electromagnetic vector current 
and charge. Solving the pertinent SD equation for the dressed quark-photon-quark (qq^y) 
vertex V 1 is a difficult and still unsolved problem, and using the realistic Ansatze for T M still 
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remains the only practical way to satisfy the WTI. The simplest solution of the vector WTI 
is the Ball-Chiu (BC) |7| vertex 



r%cW,q) = A + (q> 2 ,q 2 )^ + ^^{A^q' 2 ,q 2 )^±^ - B_{q' 2 ,q 2 )} , (3) 



where H±(q' 2 , q 2 ) = [H(q' 2 ) ±H(q 2 )], for H = A or B. This particular solution of the vector 
WTI reduces to the bare vertex in the free-field limit as must be in perturbation theory, has 
the same transformation properties under Lorentz transformations and charge conjugation 
as the bare vertex, and has no kinematic singularities. Note that it does not introduce any 
new parameters as it is completely determined by the dressed quark propagator S(q). In 
phenomenological calculations in the SD-BS approach, this minimal WTI-satisfying Ansatz 
(||) is still the most widely used qqy vertex (e.g., Refs. [[0|,^,|1^.^(J ) . A general WTI- 



satisfying vertex can be written [27] as = Tg C + AT^, where the addition Ar M does not 
contribute to the WTI, since it is transverse, (q' — q) fl AT^(q' , q) = 0. That is, Ar /i (g',g) 
entirely lies in the hyperplane spanned by the eight vectors T^q', q) transverse to the photon 
momentum k = q' — q. Curtis and Pennington (CP) advocated a transverse Ansatz for 
AY^(q',q) exclusively along the basis vector usually denoted Tg^g',^): 

AFV, q) = Tg(q', q ^f/^ ; W, q) = ^(q' 2 - q 2 ) - (q' + q)»tf - rf) . (4) 

2d(q', q) 

Then, the coefficient multiplying T$(q',q) can be suitably chosen to ensure multiplicative 
renormalizability in the context of solving fermion SD equations beyond the ladder approx- 
imation in QED 4 [|28 |. To this end, d(q',q) should be a symmetric, singularity free function 



of q' and q, with the limiting behavior lim (J 2 >>(? /2 d(q', q) = q ; for example, 

d ± (q\q) = ^^ 2 {(q' 2 ±qr+[M 2 (q' 2 ) + M 2 (q 2 )f] , (5) 

where M(q 2 ) = B(q 2 )/A(q 2 ) is the D%SB-generated dynamical mass function, which in our 
case has the large-g 2 dependence [O] in agreement with pQCD. 



The choice d = gL corresponds to the CP vertex Ansatz T^ P suggested in Ref. J28 
We will use it in analytic calculations of T(—Q 2 ,0), which are possible for Q 2 = and 
Q 2 — ► oo. However, in the numerical calculations, which are necessary for finite values of 
Q 2 7^ 0, we prefer to use the modified CP (mCP) vertex, T^p, resulting from the choice 
d — d + in Eq. (|4j). Namely, although the original CP denominator function d-(q',q) never 
vanishes strictly, the dynamical masses it contains become negligible with growing momenta, 
so that d_(q',q) can become arbitrarily small for the large spacelike quark loop momenta 
q, q' if simultaneously q' 2 — q 2 is small, causing numerical problems. On the other hand, the 



numerical calculations of T(—Q 2 , 0) employing our standard computational methods P^|JT3 



and using the mCP vertices have no such problems and they are as reliable as those using 
the minimal, BC vertices. In contrast to the BC one, the mCP vertex is consistent with 
multiplicative renormalizability by the same token as the CP one. It should also be noted 
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that r^ CP is essentially equa^j to the high-g 2 or q' 2 leading part of the vertex of Cudell et 



al. [2{J [see their Eq. (3.14), and comments below it]. 

In the present context, the important qualitative difference between the BC vertex on 
one side, and the CP as well as the modified, mCP vertex on the other side, will be that 
^bcWi <?) — > 7 M when both q' 2 , q 2 — > ±00, whereas T^ P (q J , q) — > 7^ and T^ nCP (q', q) — > 7^ as 
soon as one of the squared momenta tends to infinity. 

The correct ABJ anomaly result ([!]) cannot be obtained analytically in the chiral limit 
unless such a qqj vertex satisfying the WTI is used fT4]]Io| . We have checked by the explicit 
GIA calculation that the ir° — > 77 amplitude (|l|) is analytically reproduced employing the 
CP and mCP vertices, in the same way as when the BC vertices were employed in earlier 
applications, e.g. [fT~5|JT6|JTT|JT"^1 . 



In the case of tt°, GIA yields (e.g., see Eq. (24) in Ref. g2J) the amplitude Tff(k, k') 

T%(k, k') = -N c -i= / -^-ti{T^q — ~,k + q — Z)S(k + q - -) 
77 v ' ; 3V2J (2vr) 4 1 w 2' H T v H T 



x Y v (h + q-^q + + (* ^ k', fi ^ u). (6) 

Here, \ is the BS amplitude of both uu and dd pseudoscalar bound states: x = Xuu = Xm 
thanks to the isospin symmetry assumed here. This symmetry likewise enables us to continue 
suppressing flavor labels also on the quark propagators S and qqj vertices T M . We follow 
the conventions of Ref. |12| , including those for the flavor factors and flavor matrices A a . 
Then, Xtt°(q,p) = x(q,p) ^ 3 /v / 2, so that the prefactor 1/3a/2 in Eq. (^) is just the flavor 
trace tr(Q 2 A 3 /-\/2) where Q is the matrix of the fractional quark charges. 

In Ref. [11] (where only the BC vertex was employed), the transition form factor 
T^o^—Q 2 , 0) was numerically evaluated (for < Q 2 < 2.8 GeV 2 only) employing the soft and 
chiral limit. That is, in Ref. |TT[ we approximated (for the pion only!) the BS vertex by its 
leading, O(p ) piece: r ?r o(g,p) w T w o(q, 0) = 75 A 3 B(q 2 )/ f n . For the sake of comparison, we 
again give T 7r o(— Q 2 , 0) evaluated in that approximation (denoted by the dash-dotted line in 
Fig. 2), but now up to Q 2 m 8 GeV 2 . However, in the present work we go beyond this ap- 
proximation, using our complete solution for the BS vertex r 7r o(g,p), viz. the BS amplitude 
Xtt (q,p), given by the decomposition into 4 scalar functions multiplying independent spinor 
structures. The approximation we still keep, is discarding the second and higher derivatives 
in the momentum expansions. 

4. T 7r o(— Q 2 ,0) obtained in this way, and using the BC vertices, is depicted (after multipli- 
cation by Q 2 ) in Fig. 2 by the solid line. It tends to be on the high side of the data, but the 
general agreement of our T 7r o(— Q 2 ,0) with the experimental points is rather good, except 
for the intermediate momenta around the interval from 1 to 2 GeV 2 interval. The dashed 
line denotes Q 2 T 7r o(— Q 2 , 0) obtained in the same way, but employing the mCP vertices. For 
both curves, the agreement with experiment would be improved by lowering them somewhat 



Up to the inclusion of the mass functions M(q ), and for the simplest choice of their [29] rj- 
function: rj(q', q) = 1, i.e., their choice n = for the exponent in r/. 



(at least in the momentum region Q 2 < 4 GeV 2 ), which could be achieved by modifying the 
model |^|] and/or its parameters so that such a new solution for A(q 2 ) is somewhat lowered 
towards its asymptotic value A(q 2 — > — oo) — ► 1. (Of course, in order to be significant, this 
must not be a specialized refitting aimed only at A(q 2 ). Lowering of A(q 2 ) must be a result 
of a broad fit to many meson properties, comparable to the original fit |[22|| . This, however, 
is beyond the scope of this letter.) 

While the empirically successful anomaly result ([!]) at Q 2 = is model-independent 
in any consistently coupled SD-BS approach (being the consequence only of the correct 
chiral behavior due to the incorporation of D^SB), the transition form factor for Q 2 ^ of 
course depends not only on the WTI-preserving gg7-vertex Ansatz, but also on the chosen 
bound-state model. However, it must be stressed that we we did not do any parameter 
fitting whatsoever, as we have used the parameters obtained from Jain and Munczek's | 22[ 



broad fit to the meson spectrum (except the rj-rj' complex) and the pseudoscalar decay 
constants. Our model choice [^2j has subsequently given us predictions for r] c , rjb — > 77 
[TT| , the description of rj-rf complex and 77, rj' —>■ 77 [[12] , and now, with the same set of 
parameters, the transition form factors which are empirically successful for the presently 
accessible high values of Q 2 , and at the same time in agreement (in the chiral limit even 
exactly and analytically |T^-|rB|,|TT],|r2[ ) with the anomalous amplitude (|XJ) at Q 2 = 0. 

Our transition form factors thus also agree rather well with other successful theoretical 
approaches. The vector-meson dominance (VMD) model (which is still the most successful 
one from the purely phenomenological point of view ||) and the QCD sum rule approach 
H give the transition form factor which is some 10% below our "BC" T 7r o(— Q 2 , 0) for large 
values of Q 2 . Therefore, in the large Q 2 region, our "BC" values are halfway between the 
uppermost line in Fig. 2 denoting the asymptotic pQCD Jl8| version {J = 2) of the result 
(0), and VMD as well as the QCD sum rule results ]|] of Radyushkin and Ruskov [amounting 
to Eq. (HD with J « 1.6]. By the highest presently accessible momenta, our "mCP" version 
of T 7r o(— Q 2 , 0) has crossed to the low side of both other theoretical predictions and the data, 
but is still within the error bars. 

The large-Q 2 leading power-law behavior (121) was first derived from the parton pic- 



ture in the infinite momentum frame - e.g., see Ref. ||18|| . In this and other similar 
pQCD approaches, the precise value of the coefficient of the leading l/Q 2 term depends 
on the pion distribution amplitude (p % (x) which should contain the necessary nonpertur- 
bative information about the probability that a partonic quark carries the fraction x of 
the total longitudinal momentum. E.g., the well-known example of a "broad" distribution 
ip% z (x) = f n 5v3x(l — x)(l — 2x) 2 (proposed by Chernyak and Zhitnitsky motivated 
by sum-rule considerations) leads to J = 10/3, but it is too large in the light of the latest 
CLEO data |]. In contrast, the asymptotic (f^(x) = f n y/3x(l — x) favored by Lepage and 
Brodsky [18] yields J = 2, resulting in the line touching the error bars of the presently high- 
est Q 2 data from above in Fig. 2. In fact, in the strict InQ 2 — > 00 limit, every distribution 
amplitude must evolve into the asymptotic one, <p^{x) — > f^(x), if the effects of the pQCD 
evolution are taken into account. However, even at Q 2 - values considerably larger than the 
presently accessible ones, other effects may still be more important than the effects of the 
pQCD evolution. This is the reason why other approaches and other forms of (p n (x) should 
be considered even when they do not incorporate the pQCD evolution. The form 
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(7) 



is suitable for representing various distribution amplitudes because it is relatively general 
fl~8|| : ( > 1 yields the (empirically favored) distributions that are "peaked" or "narrowed" 
with respect to the asymptotic one (( = 1), whereas ( < 1 gives the "broadened" distri- 
butions which, however, now seem to be ruled out for the same reason as (p^ z (x) quoted 
above, since J > 2 is ruled out empirically by CLEO IJ. Namely, it is easy to see that 
Eq. (0) implies ( = 2/(3J-4). 

Indeed, at the highest presently accessible momenta, the asymptotic prediction (J = 2) 
is lowered by some 20% by the lowest order QCD radiative corrections |JIJ, amounting to 
J pa 1.62, which fits the CLEO data well. Of course, these 0(a s ) corrections mean that 
Q 2 T 7T o(—Q 2 ,0) is not strictly constant, but according to Ref. ffl] it rises towards 2f w so 
slowly (by just 4% from Q 2 = 9 GeV 2 to Q 2 = 36 GeV 2 ) that we can take it constant in 
practice. The situation is similar with the sum-rule approach of Radyushkin and Ruskov 
0. Their Q 2 T n o(— Q 2 , 0) starts actually falling after Q 2 ~ 7 GeV 2 , but so slowly that Eq. 
(§) with the constant J pa 1.6 represents it accurately at the presently accessible values of 
Q 2 . This corresponds to a rather narrow distribution (|7]) with ( = 2.5. 

The same type of the leading large-Q 2 behavior as in Eq. (H), was obtained by Manohar 
| 20| using the operator product expansion (OPE). According to his OPE calculation, the 
coefficient in Eq. (0) giving the leading term is J = 4/3, which is below our large-Q 2 
T 7r o(— Q 2 ,0) by ~ 20% when we use the BC vertex, but - as we will see below - exactly 
coincides with the Q 2 — > oo limit obtained using the mCP and CP vertices. The coefficient 
J = 4/3 is the lowest one still consistent with the form (0) because it corresponds to 
( = oo. The pion distribution amplitude (|?P then becomes infinitely peaked delta function: 

¥v(x) = (/ w /2v^)5(x-l/2). 

For Q 2 > 4 GeV 2 , our "BC" T 7r o(— Q 2 ,0) also behaves in excellent approximation as 

Eq. (||), with J p^ 1.78. This J would, in the pQCD factorization approach, correspond 

to (f n (x) (|7|) with ( = 1.5. On the other hand, our "mCP" T 7r o(— Q 2 ,0) falls off faster than 

1/Q 2 for even the largest of the Q 2 values depicted in Fig. 2. However, it does not fall off 

much faster, as our "mCP" Q 2 T 7r o(— Q 2 , 0) at Q 2 = 18 GeV 2 is only 6%, and at the huge 

Q 2 = 36 GeV 2 is only 10 % smaller than at Q 2 = 9 GeV 2 (i.e., at the highest presently 

accessible momenta). Moreover, we can show analytically that, generally, the l/Q 2 -behavior 

of Eq. (0) must at some point be reached in our approach, although Fig. 2 shows that for 

the mCP qq'j vertices this can happen only at significantly higher Q 2 than it happens for 

the BC qqj vertices. 

5. It is very pleasing that in the present approach (the coupled SD-BS approach in conjunc- 
tion with GIA), besides numerical results, one can get also some analytical insights in the 
regime of asymptotically large Q 2 . This way we can make more illuminative comparisons 
with the asymptotic behaviors predicted by the various approaches - notably pQCD on the 
light cone |1§[19| and OPE @]. 

Since the pion is light, k ■ k' pa Q 2 /2 for large negative k 2 = —Q 2 and k' 2 = 0. Taking 
into account the behavior of the propagator functions A(q 2 ) and B(q 2 ) (e.g., see |13[]) we can 
then in Eq. @ approximate those quark propagators that depend on the photon momenta 
k and k' , by their asymptotic forms: 
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%-! + *)« - 1) « ~Qt - ijO , (8) 

and analogously for the propagator where k is replaced by k' . Although the relative loop 
momentum q can be large in the course of integration, its neglecting is justified because the 
BS amplitude xijiiP) decays quickly and thus strongly damps the integrand for the large q's. 

For the moment let us make an additional approximation by replacing the dressed elec- 
tromagnetic vertices in Eq. (|6|) by the bare, free ones: r^(g, q') — * 7^. The tensor amplitude 
(|) then becomes 

T^(k, k') = ^J yj^tr ( x (g,p) [(* - |) a7 ^V + ~ f )a 7 ^V]) , (9) 
which is readily rewritten as 

T%(k,k') = --pe^(k-k') x j -0^tr[ X (q,Ph,l 5 ] = e^(k - k') xPa U ■ 

(10) 

The second equality holds since in the Mandelstam formalism, the integral in Eq. (|T0|) is 
equal to — ip a f n \[2jN c - e.g., see |2l|Jl2]1 . Using the definition of the scalar amplitude 



T 7r o(/c 2 , k' 2 ), we finally recover Eq. (|2|) with the coefficient having the special value J = 4/3. 
We have thus found that the asymptotic behavior predicted by the present approach with 
bare qq^y vertices (and, as seen below, also with dressed mCP and CP vertices), is in exact 



agreement with the leading term predicted by OPE |20| . 

The asymptotic behavior oc l/Q 2 obtained in the present approach is especially satisfying 
when compared with that resulting from the calculation of the triangular quark loop carried 
out in the simple constituent quark model (with the constant light-quark mass parameter 
m u ), where T 7r o(— Q 2 ,0) oc (rn 2 L /Q 2 )\n 2 (Q 2 /ml) as Q 2 — ► 00, which overshoots the data 
considerably || because of the additional ln 2 (Q 2 )-dependence. 

As seen from Fig. 2, the asymptotic behavior Q 2 T n o(— Q 2 , 0) = (4/3)/^ barely touches 
the experimental error bars from below. However, Manohar |20| pointed out that his OPE 
approach also indicates the existence of potentially large corrections to his leading term. 
(Note that he also pointed out that, for the same reason, significant corrections should 
similarly affect the light-cone prediction {J = 2) of Lepage and Brodsky |T8]] .) It is easy 
to see that in our case when the BC vertex is used, the main origin of the corrections 
raising J from J = 4/3 resulting from the usage of the bare vertices, to J ~ 1.78 of 
our full numerical calculation, is the usage of these dressed qqj vertices. It is nevertheless 
instructive to formulate another successive approximation to illustrate the gradual transition 
to the T n o(—Q 2 , 0) resulting from the full calculation. The idea is that at high Q 2 , only the 
iJ + -term contributes significantly in the BC- vertex (|j), as the if_-terms are suppressed by 
their Q 2 -dominated denominators in the both BC- vertices: pf — p% = ±(g =(= p/2) 2 =F (k + 
q — p/2) 2 w ±Q 2 /2. The detailed derivation shows that, as Q 2 grows, this indeed happens, 
and only ^^[A^pD + A(pl)} contributes significantly to both of the BC vertices. Since p\ 
in one of the vertices, and p\ in the other, also contain Q 2 /2, the BC vertices reduce to 
7 M [1 + A([q ±p/2} 2 )]/2 since A(-Q 2 /2) -»• 1 as Q 2 -> 00. The adequacy of the chiral limit 
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approximation p 2 = M 2 = (and, in the case of our particular model choice, of neglecting 
the derivatives of A in the expansion of A([q±p/2} 2 ) around q 2 since our solution for A(q 2 ) 
is a smooth, slowly varying function ||13|| ) leads to the asymptotic expression 

3Q 2 



T.o(-Q 2 ,0)^-^ (11) 



where f n is given by the same Mandelstam- formalism expression as f n , except that the 
integrand is modified by the factor [1 + A(q 2 )] 2 /4. A further approximation would be to 
substitute A(q 2 ) by some average value, say A((f) where q 2 is a value from the — g 2 -interval 
[0.5, 2] GeV 2 of the momenta that have the most influence on the integral that gives f n . The 
relation f n = f n [l + A(g 2 )] 2 /4 can serve as a definition of q 2 . The asymptotic value can thus 
- in the case of the BC vertex (H) - be expressed as Eq. @ with J — [1 + v4(g 2 )] 2 /3. In the 
case of our solutions, the slow and moderate variation of the shape of A(q 2 ) (see Fig. 3 in Ref. 
|j~3|l ) permits the approximation A((f) « (A min + A max ) / 2 = (l+A max )/2. In our special case 
this can be further approximated by A(0), leading to the estimate Alq 2 ) ~ ^4(0) = 1.25 |13 



This approximation would imply J »s 1.69, whereas the more accurate Eq. ([11]) through 
f n ~ 1.334/^ leads to Eq. (|2|) with J « 1.779, which is in fact numerically indistinguishable 
from the asymptotics of our full calculation and the empirical fit. 

However, all this happens for the BC vertex, where the "soft" leg adjacent to the pion 
BS-amplitude always contributes A(q 2 ), but not for the mCP vertex, nor for the CP one. 
Since they tend to the bare vertex as soon as the high momentum flows through one of its 
legs, f-x does not get replaced by f n , and the "bare" result, Eq. (^) with J = 4/3, continues 
to hold for Q 2 — > oo when the mCP or CP vertex is used. Of course, the usage of the mCP 
vertex or the CP vertex instead of 7^ causes considerable differences for the finite Q 2 . 
6. The most important result of the present paper is that we have demonstrated that 
the modern version of the constituent quark model which is given by the coupled SD-BS 
approach, provides from Q 2 = to Q 2 —>■ 00 the description for 77* — > 7r° which is - 
independently of model details - consistent with the Abelian ABJ anomaly and with the 
QCD predictions (0) for the leading large-Q 2 behavior. Since f n is a calculable quantity 
in the SD-BS approach, the model dependence is present in the (successfully reproduced 
P^ , PT| -P^|) value of f n , but our derivation of the asymptotic forms (0) and flTT| ) is model- 
independent. Of special importance is also that this derivation [Eqs. (|9|), ([It])] of the 
large-Q 2 behavior applies to both Minkowski and Euclidean space. The same holds for the 
considerations involving /„. and the [1 + A(g 2 )] 2 /4 correction factor, which have enabled us 
to make smooth and accurate contact between our large-Q 2 numerical results calculated 
with the Euclidean solutions of the chosen model |[22|| , and analytical results on the large- 
Q 2 asymptotics. The agreement between the analytical and numerical results enhances 
our confidence in the accuracy of our numerical methods and procedures (employing the 
Euclidean bound-state solutions) in this and earlier papers |~1]-[13| • 



To mention one possible application, the relevance of the present work for the light-by- 
light scattering contribution to the muon g — 2 is obvious. Nevertheless, the relevance of its 
extension to the case when both photons are off-shell, i.e., to the 7r°7*7* form factor, is even 
larger (for g — 2), since there are no experimental data on T 7r o(— Q 2 , —Q' 2 ) at present. In 
the light of its consistency with OPE and axial anomaly, but also the (at least) qualitative 
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consistency with the QCD sum rules, light-cone pQCD and VMD dominance, we propose 
that the present approach be used in a re-calculation of some 7r°7*7*-dependent predictions 
of Ref. || instead of the double- VMD Ansatz (or instead of the Extended Nambu-Jona- 
Lasinio model which does not give the correct asymptotic behavior even for the 7r°7*7 
vertex |3l|). The reason is that this VMD- motivated Ansatz for the 7r°7*7* form factor 
(Eqs. (3.5) and (4.1) in Ref. 0) disagrees with what it should be according to the simple 
and unambiguous extension of T n o(— Q 2 ,0) presented in this work. In particular, the VMD 
Ansatz for T 7r o(— Q 2 , — Q' 2 ) behaves as 1/(Q 2 Q' 2 ) for asymptotically large Q 2 ,Q' 2 , whereas 
it should be only as the inverse of the sum of these squared momenta. Our full treatment 
of the 7r°7*7* form factor will be given in another, more detailed paper ||32|| , but already 



here we can point out that it is almost trivial to generalize the derivation of the asymptotic 
behaviors (EJ) and ( 11 ) to the case when both photons are off-shell, k 2 = —Q 2 << and 



k' 2 = —Q' 2 < 0. Not much changes, except that k ■ k' {Q 2 + Q' 2 )/2. Our prediction with 
the bare (or mCP, or CP) (7(77- vert ices for this generalization of Eq. (EJ) turns out to be 



TA-Q\-Q*) = l^, (12) 



which is in agreement with the leading term of the OPE result derived by Novikov et al. 
33fl for the special case Q 2 = Q' 2 . The distribution-amplitude-dependence of the pQCD 



factorization approach cancels out for that symmetric case, so that T n o{—Q 2 , —Q' 2 ) in this 
approach {e.g., see p4||), in the limit Q 2 = Q' 2 — > oo, exactly agrees with both our Eq. 



(|12|) and Ref. The inclusion of the dressed, WTI-preserving vertices does not require 
any modifications of the asymptotic forms in the case of qq'j vertices (such as mCP and CP 
ones) which reduce to the bare one as soon as one of their quark momenta squared tends 
to infinity. In the case of the BC vertex (|3|), the generalization of the better-approximated 
expression (|TTf) requires nothing except the substitution f n — > f n in Eq. fll2|) . For our [13| 



particular A{q 2 ), f n can, in a good approximation, be factored as explained above. 

We note that the VMD-motivated Ansatz for the P°7*7* {P° = 7r°, 77, 77') transition 
form factor is also used as a theoretical input in the CLEO data analysis - see Eq. (8) 
in Ref. |[J]. A re-analysis of the data in the spirit of the present insights from the SD- 
BS approach, is therefore desirable. At least, one should examine the effects on the data 
analysis when one makes a change in Eq. (8) of Ref. § like 1/[(1 + Q 2 /A 2 ){1 + Q' 2 '/ 'A 2 )} -> 
1/(1 + Q 2 /A 2 + Q'7A 2 ). 

Actually, 77 and rf appear not only in the CLEO data ||, but are also relevant for the 
g — 2 analysis || discussed above. This underscores the importance of extending the present 
work to include the 777*7 an d 7/7*7 transition form factors. While the full calculation thereof 
obviously must be relegated to another paper |32| , our predictions for the Q 2 —>■ 00 behavior 



of the SU(3)/ states 77 8 and 770, and also 77 c and r/b, can be easily obtained by redoing the 
derivation of the asymptotic expression so that the track of their flavor content is kept. But, 
one then realizes that the generalization to both photons being virtual can be easily done 
also for r)8,r)o,r] c and r] b in the way discussed two passages above for 7r°. Therefore, in the 
approximation of the bare qq^ vertices, in the asymptotic regime with one of the photon 
momenta large (say Q 2 — > 00), and the other (say Q' 2 ) anywhere between and 00, the 7/g 
and 770 transition form factors are given by the expressions analogous to Eq. ([12]) , except 
that its factor 4/„-/3 = C w o, appropriate for P° = 7r°, is replaced by the respective factors 
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C m = (4/9\/3)(5f n - 2f aS ) and C Vo = (4 v / 2/9 v / 3)(5/ 7r + f ag ). The auxiliary quantity f sS is 



the decay constant of the unphysical ss pseudoscalar bound state, evaluated in [[12] to be 
f S s = 1.47/tt for the chosen bound-state model [Q, yielding the model-dependent values 
C m = 0.53 U = 49.3 MeV and C no = 2.35 f v = 219 MeV. 

When the r/s^o asymptotic expressions are improved by using the dressed, WTI- 
preserving vertices instead of the bare vertex 7 M , the asymptotic expressions do not change 
at all in the case of the mCP and CP vertex, and in the case of the BC-vertex the only 
changes are again just the substitutions f n — > f w and / s5 — > f sS - 

Similarly, for r] c and rjb, the analogous expressions when at least one of Q 2 or Q' 2 is 
large, are obtained from Eq. (|T2D by replacing 4/ 7r /3 = C^o with C Vc = (16v^2/9)/r, c and 
C Vb = (4^/9) f Vb , respectively. (Of course, if Q' 2 < Mp , neglecting the meson mass M P o 
is for Q' 2 < M| a rougher approximation than setting Q' =0. M Vc and especially M Vb are 
already so large that keeping k ■ k' = (Q 2 + Q' 2 + Mj, )/2 in Eq. (|T^) can be important for 
momenta Q 2 accessible in practice.) In our chosen bound-state model P2| , we have obtained 
f Vc = 213 MeV and / % = 284 MeV, leading to the model-dependent values C Vc = 536 MeV 
and C Vb = 179 MeV. This is of immediate importance for experiment, as the r] c transition 
form factor T Vc (— Q 2 ,0) can be measured at L3 |35]. 
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FIGURE CAPTIONS 



Fig. 1: The diagram for P° — > 77 decays (P° = tt°, 77, rj', rj c , rjb). Within the scheme of gener- 
alized impulse approximation, the propagators and vertices are dressed. 

Fig. 2: The comparison of our results for Q 2 T n o(— Q 2 , 0), the pion transition form fac- 
tor, with the CELLO (circles) and CLEO (triangles) data and with the Brodsky- 
Lepage interpolation formula (the dotted line, which also denotes the pQCD limit- 
ing value 2f n ). The dash-dotted line represents our Q 2 T 7r o(— Q 2 , 0) evaluated as in 
our earlier paper [11] (but now to higher squared momenta), i.e., exclusively with 
the BC qq^ vertices, and in the chiral and soft limit approximation for the pion 
\T^{q,p) « r^o(g,0) = 7 5 A 3 J B(g 2 )//^]. Our results for Q 2 T^{-Q % , 0) without that 
approximation are depicted by the solid line for the case of the BC vertices, and by 
the dashed line for the mCP vertices. The one obtained with the BC vertex practically 
saturates beyond Q 2 ~ 6 GeV 2 at higher values [Q 2 T n o(— Q 2 , 0) w 164 MeV] than that 
obtained with the mCP vertices. This latter one gets lower beyond Q 2 ~ 2.5 GeV 2 
and does not yet saturate at the presently accessible momenta although approaches 
asymptotically Q 2 T n o(—Q 2 , 0) — > 4:f n /3 for much higher Q 2 . (This limit, which is also 
the leading term resulting from OPE [p(| , is denoted by the dashed straight line.) 
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